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Abstract. - When a liquid drops impinges a hydrophobic rough surface it can either bounce 
off the surface (fakir droplets) or be impaled and strongly stuck on it (Wenzel droplets). The 
analysis of drop impact and quasi static "loading" experiments on model microfabricated sur- 
faces allows to clearly identify the forces hindering the impalement transitions. A simple 
semi- quantitative model is proposed to account for the observed relation between the surface 
topography and the robustness of fakir non-wetting states. Motivated by potential applica- 
tions in microfluidics and in the fabrication of self cleaning surfaces, we finally propose some 
guidelines to design robust super hydrophobic surfaces. 



Some plants leaves and insects shells exhibit extreme hydrophobicity, making the deposi- 
tion of water drops on their surface almost impossible [1]. All these superhydrophobic biosur- 
faces share two common features: they are made of (or covered by) hydrophobic materials, 
and are structured at the micron and sub-micron scales. 

During the last decade much effort has been devoted to design artificial solid surfaces with 
comparable water-repellent properties. Their potential applications range from lab on a chip 
devices to self cleaning coating for clothes, glasses,... The actual strategy consist in mimicking 
superhydrophobic biosurfaces designing rough substrates out of hydrophobic materials. To 
achieve this goal both top-down and bottom-up approaches have been successfully developed: 
chemical synthesis of fractal surfaces [2], growth of carbon nanotube forests [3], deep silicon 
dry etching [4], see also [5] and references therein. We briefly recall the paradigm to account 
for superhydrophobicity. Two different wetting states can be observed on microstructured 
hydrophobic surfaces: (i) Wenzel state: the liquid follows the topography of the solid surface. 
Defining the surface roughness ( as the ratio between the total surface area over the apparent 
surface area, the equilibrium contact angle of a liquid drop is given by cos((9) = C cos (#flat) 7 
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Fig. 1 - Rapid camera snapshots: (a) Impact of a water droplet on a micropatterned surface: Pattern 
geometry Si(h = 8/im), impact speed V\ — 0.27 m.s -1 , droplet radius 1mm. The droplet bounces 
off the surface, (b) Impact on the same substrate for V\ — 0.6 m.s -1 , droplet radius 1 mm. The water 
droplet remains stuck on the surface. Time intervals ^15 ms for the two sequences. 

where #fl a t is the Young contact angle on the flat surface [7]. (ii) fakir state: The liquid only 
contacts the highest parts of the rough solid, air pockets remain trapped between the solid 
and the liquid surface. Only a fraction <p of the solid surface, corresponding to the extremities 
of the protrusions, is wetted by the liquid. As proposed by Cassie, a water drop adopts 
a contact angle given by the weighted sum cos#f a ki r = </>cos#fl at — (1 — 0), (Cassie-Baxter 
relation) [8] . More precisely, it has been shown that the surface energy is minimal in the fakir 
(resp. Wenzel) regime if f is larger (resp. smaller) than cos#f a ki r / cos#fl at . In other words, 
the rougher the substrate, the more the fakir state is energetically favored. 
The characterization of superhydrophobic surfaces is usually restricted to equilibrium contact 
angle measurements. However, it has been recently reported that: (i) the value of the measured 
contact angle strongly depends on the way the droplet is deposited on the surface [9, 10] 
and (ii) droplets squeezed between two moderately rough surfaces can undergo a sharp and 
irreversible transition from a fakir to a lower contact angle Wenzel [9]. Moreover, the drops 
do not only reduce their contact angle, but also increases their contact angle hysteresis, the 
contact line appears to be strongly pinned on the substrate, any self-cleaning properties is 
thus definitely lost. Today, a clear picture allowing to identify the mechanisms responsible 
for the impalement transition from a fakir to a Wenzel wetting state is missing. Only few 
theoretical explanations have been attempted [12-14], see also [15] in a different context. 
Beyond the usual thermodynamic approach, it appears crucial to extend the understanding of 
such impalement transition to drop impact dynamics. Indeed almost all practical applications 
of superhydrophobic surfaces rely ultimately on their ability to repel impinging drops (rain 
drops, sprays,...) 

In this Letter, we present the characterization of the drop impalement transition on micro- 
fabricated surfaces under quasi-static and impact dynamics as well. The dual analysis of the 
two series of experiments is completed by a simple model that allows to unambiguously iden- 
tify the forces hindering drop impalement and to propose a unified criterion for the robustness 
of fakir non- wetting states. 

All the presented experiments have been performed on PDMS silicon elastomer (poly- 
dimethylsiloxane, Sylgard 184 Dow- Corning) surfaces, micropatterned using classical soft- 
lithography molding methods [11]. The microfabricated surfaces are triangular arrays (pitch 
p) of cylindrical pillars (radius r). Varying the thickness of the primary mold made of pho- 
toresist resin (SU08, Michrochem), we have varied the pillar height from 2.7 /im up to 75 /im. 
Two different patterns have been used: Si (r = 11 /im, p = 50 /im) and 62 (r = 9 /im, 
p = 42 /im). Both patterns have the same pillars density: <j) ~ 0.15. Water drops lying on flat 
PDMS surfaces have an advancing (resp. receding) contact angle # a = 110° (resp. T = 80°). 
When gently deposit on the patterned surfaces, the measured contact angles of the water 
drops agree with the Cassie-Baxter relation (advancing angle a ~ 155°). Fakir states are 




D. BARTOLO et all BOUNCING OR STICKY DROPLETS: IMPALEMENT TRANSITIONS 



3 



N0N " BOUNCING 
BOUNCING 



(a) 



0.4 0.8 
Impact velocity (m.s" 1 ) 





I I I > 

I, I 


STICKY 


y 


BOUNCING 












NON-BOUNCING " 



40 
h (urn) 



Fig. 2 - (a) Inverse of the impact contact time plotted versus the impact velocity, V\. Surface 
geometry: pattern S±, pillars hight, h = 9/xm. (b) Critical impalement kinetic energy plotted versus 
pillars height, for Si pattern geometry. A: Critical energies delimiting the Bouncing-Sticky transition. 
Plain line: linear best fit of the data (0 < h < 20 /xm. Dashed line: best constant fit (h > 20 /mi), y: 
critical energies delimiting the bouncing-non-bouncing transition. Dotted line: best constant fit. 



hence observed despite their surface energy would be minimized in the impregnated Wenzel 
state (except for h = 75 /im pillars that should insure equilibrium in the fakir state). 

A first experiment consist in studying the impact of water droplets (radius 1 mm) delivered 
by a precision needle on the micropatterned PDMS substrates. Increasing the fall height 
increases the impact velocity V\. The impact events have been observed using a high-speed 
video system (frame rate: 1000 fps). 

The impinging drops first expand rapidly. Subsequently, due to the hydrophobicity of 
the surface, the drops retract and sometimes bounce off the surface as illustrated in Fig. [I] 
In Fig. 0a., the inverse of the contact time is plotted versus the impact speed. We observe 
that bouncing occurs only in a range of impact velocities (Vnb < V\ < Vbs)- As previously 
discussed by Richard et at. in [16], the contact time does not depend on V\. Looking more 
carefully at the retraction dynamics, three distinct regime can be identified. 

- V\ > Vqs : sticky droplets. In this regime, the contact line hardly retracts and the 
intantaneous contact angle reaches values as small as 40°, see Fig^b. This strong pinning is 
a clear evidence that the microstructure has impaled the liquid surface. Observations with a 
microscope have systematically confirmed that the pillars are impregnated. 

- Vnbb < Vj < Vbs' bouncing droplets. At intermediate velocities, the drop bounces on 
the surface (this behavior is never observed on flat PDMS), several bouncing events can be 
observed. Finally, the drop remains on the surface adopting a large contact angle consistent 
with the Cassie-Baxter prediction. The initial kinetic energy of the drop is not sufficient to 
overcome the energy barrier hindering the impalement transition. 

- Vj < Vnbb : non-bouncing droplets. As previously reported in [17], we obseve a low 
speed threshold below which droplets do not bounce anymore. In this regime, the drop 
weakly expand after impact. Then, the drop undergoes damped oscillations to reach its a 
quasi-spherical shape, corresponding again to a fakir non-wetting state. Though bouncing is 
not observed, water does not fill the microstructure. 

Varying the height of the pillars we construct the phase diagram plotted in Fig. Ob. The 
kinetic energy thresholds delimiting the three regimes are plotted versus h. 

The non-bouncing to bouncing transition can be easily understood. During the drop 
retraction stage, we assume that the pinning of the contact line is the main source of kinetic 
energy dissipation. Denoting 7 the liquid-air surface tension and Acos# the contact angle 
hysteresis of the fakir drop (independent of h), we can assess the pinning force per unit length: 
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Fig. 3 - (a) Two consecutive pictures of a compressed droplet viewed from above. The light part inside 
the droplet is in contact with the microstructured substrate. Between the two snapshots the droplet 
has transited from a fakir to a Wenzel state. Note that the contact line on the Si surface adopts 
the hexagonal symmetry of the pattern. External drop radius (b) Three consecutive pictures taken 30 
sees apart of an evaporating droplet gently deposited on a substrate Si, (h = 20 /mi). Between the 
two last snapshots the droplet has transited from a fakir to a Wenzel state. Contact radius of the 
drop on the last picture:450 /im 



7A cos 6. The minimal kinetic energy needed to observe bouncing thus scales 7P 2 A cos 6. The 
order of magnitude of this energy threshold is ~ 5 10 _T N.m, in good agreement with what is 
experimentally observed. 

The critical impalement threshold delimiting the bouncing and the sticky regimes exhibits 
nontrivial variations with h, see Fig. 0b. For short pillars, the critical kinetic energy increases 
linearly with h. Above h = 20 /im, it becomes independent of the texture roughness. We insist 
on the irreversible nature of the impalement transition whatever the equilibrium wetting 
states. This result is one of our most important finding. We now have to identify the physical 
mechanism hindering the liquid impalement upon impact of a droplet. A priori, both capillary 
and hydrodynamic forces, act on the liquid surface to impede the impalement process. We 
disentangle the two effects by performing two sets of systematic quasistatic experiments. This 
allows to asses to role of the sole capillary forces. 

- Squeezing [18]: A droplet is squeezed between, a microstructured substrate and a fluo- 
rinated glass slide (advancing contact angle with water 120°). The drop is observed through 
the glass plate (see Fig. |3Ja). The gap between the substrate and the glass plate is slowly 
decreased until we observe a rapid jump forward of the contact line on the PDMS substrate. 
This jump is the signature of the liquid impalement. The increase of the drop Laplace pressure 
is the motor behind the liquid impalement. Fitting the droplet shape just before the transi- 
tion by a surface with a constant mean curvature (C), we infer a critical impalement Laplace 
pressure Pi mp = 7C. Fig. 0] For the microstructures made of long pillars, our experimental 
setup did not allow to reduce the gap between the two solid surfaces sufficiently to observe the 
impalement transition. To circumvent this technical obstacle, we used an alternative method 
to increase the Laplace pressure in the drop. 

- Evaporation: A millimetric water droplet is gently deposited on the microstructured 
substrates. Since the drop evaporates, its radius slowly decreases with time, hence its curva- 
ture and the Laplace pressure pushing the liquid surface on the micropillars rise continuously. 
After few minutes, the drops adopts its receding contact angle (~ 120°) and the contact line 
retracts. Again, above a critical pressure Pimp, one can observe a sudden variation of the 
contact angle and a strong pinning of the contact line that stops its retraction until complete 
evaporation, see Fig. |31b. These two observations witness the drop impalement transition. 
Once a drop has reached a Wenzel state, it never relaxes toward a fakir conformation. We 
emphasize that arbitrarily small drops cannot be maintained in a fakir state. This results 
seems to be at odd with the the criterion: £ > cos#f a ki r / cos#fl at which is independent of 
the drop size. We point out that this criterion has been established ignoring any pressure 
difference across the liquid interface, i.e. for infinitely large drops. 
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Fig. 4 - Impalement pressure Pimp plotted versus the pillars' height h. Circles: evaporation experi- 
ments, squares: compression experiments, triangles: drop impact experiments. Pattern geometry: Si. 
Plain line: best linear fit (h < 20 /xm). Dashed line: best constant fit (/i > 20 /mi). Inset: Impalement 
pressure normalized by -Pcrit* mg > see Eq. |3] plotted versus the pillars' aspect ratio h/r. Filled (resp. 
open) symbols correspond to surface Si (resp. S2). 



To quantitatively compare the outcome of our quasistatic and impact experiments, we 
assess the pressure pushing the liquid interface when a water drop hits the micropillars. 
Neglecting any viscous effect, the dynamic pressure acting on the liquid interface scales as 
P ~ \pV\ \ with p the liquid density. Fig. 21 gathers static and impact experimental data. 
The collapse of our data on a single master curve is a strong evidence that, hydrodynamic 
forces do not play any significant role in the impalement transition on our surfaces rough at 
the 10 microns scale. The wetting of the microstructure is mainly hindered by capillary forces. 
This constitutes our second main results. 

We have unambiguously identified the forces impeding the fakir to Wenzel transition, we 
can now propose a simple semi-quantitative model to account for the relation between the 
pattern geometry and the critical impalement pressure. Whatever its precise origin (curvature, 
liquid flow,...), the internal drop pressure pushes the liquid interface downward. A force 
Fp = P [A(l — (j))] is experienced by the liquid free surface enclosed by an elementary cell of 
the lattice, where we have denoted A the projected area of one cell. At mechanical equilibrium, 
this force is balanced by the capillary force, Fq = N p [27rjr cos(0)] applied at the top of 
N p = (j)A/7rr 2 pillars, see Fig.0 Here is the "average" contact angle defined on the pillars 
sides. Writing explicitly the equilibrium condition Fp + Fq = 0, we obtain: 

P=-^-|cos(0)| 2 . (1) 
1 — r 

Note that in the above equation the precise shape adopted by the liquid interface is encoded 
solely in the cos(#) prefactor. 

We now propose two impalement scenarios characterized by two critical values of the 
contact angle on the pillars. 

- "Touch down" scenario: increasing the drop pressure, the curvature of the interface 
increases. This implies that the minimal height separating the liquid interface and the basal 
surface of the substrate diminishes. Fakir states cannot be stabilized if this minimal height 
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Fig. 5 - (a) Sketch of the "touch down" scenario for the fakir-sticky droplet transition, (b) Sketch of 
the "sliding" scenario for the fakir-sticky droplet transition 



goes to zero, see Fig. 5. a. This contact condition can be expressed in term of a critical angle 
value 0(h). An exact computation of this critical angle would require the full determination of 
the drop shape. To bypass this difficulty, we can estimate this critical angle in a much simpler 
manner. Using a the small deformation approximation it is straightforward to compute profile 
of a fluid interface lying on top of two concentric cylinders (radius r and p) . For this simplified 
geometry, one can easily show that, cos[#(ft)] ~ h/r, for r/p <C 1 and omitting logarithmic 
corrections. It then follows that the critical "touch down" impalement pressure scales at first 
order with respect to the solid fraction <\> [19]: 



imp 



7ft 
P Z 



(2) 



where we have used (j) = (2tt / \^) (r / p) 2 for a triangular lattice. Note that for denser patterns, 
i.e. for 2r c± p, cos [0 (ft)] woulf scale as ~ ft /p. The above equation would then be modified 
and i^p tact ~ jhr/p 3 . For weakly rough substrates, Eq.0 correctly predicts the linear scaling 
of the impalement pressure with respect to the pillars height reported in Fig. and Fig. 0] 
Though several simplifications were made one can compare the estimate of the numerical 
prefactor to its experimental value. For the Si surface, Eq. predicts a slope p contact /^ ~ 
30 Pa./im -1 , the slope extracted from Fig. 0]is 2 times smaller. We do not think that this slight 
difference arises solely from the approximations needed to estimate 0(h), a precise numerical 
description of the free surface actually increases this discrepancy. We will show elsewhere that 
localized heterogeneities are good candidates to account for the overestimate of the p contact /ft 
value [20]. Indeed, the effect of chemical or geometrical heterogeneities cannot be captured 
by our simple "mean field" model. 

- "sliding" scenario: For higher pillars a transition to a critical impalement pressure inde- 
pendent of the aspect ratio h/r can understood as follow. The contact angle has another 
upper limit, namely the local advancing contact angle value a . If exceeds a , the contact line 
will spontaneously slide downward along the pillars to reach the floor of the microstructure. 
Therefore, this critical "sliding" impalement pressure is obtained taking = a in Eq. 



-^-\cos(e a )\l. 

1 — r 



(3) 



Beyond the correct scaling prediction equation Eq. |3| provides a rather good estimate of the 
experimental prefactor value for the two tested patterns, see inset in Fig. EI A factor of 2 is 
as good as we could have expected given the simplicity of our model and the precision of the 
microfabrication process. 
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This simplified model conveys a clear picture to account for the relation between the mi- 
cropattern geometry and the robustness against liquid impalement. The drop will undergo 
an impalement transition if the pressure in the drop exceeds min(i^^ p , i^p)? the crossover 
between the two impalement scenarios corresponds to an aspect ratio h/r ~ 1 for low 
patterns. 

We eventually emphasize that the above description may be a useful guide for the de- 
sign of super hydrophobic surfaces. Efficient water repellent surfaces must obviously exhibit 
a high Young contact angle, a small contact angle hysteresis and a strong resistance against 
irreversible impalement. The two first requirements can be achieved reducing the solid frac- 
tion (j). The sole comparison of the surface energies associated with Fakir and Wenzel states 
would lead to the fabrication of substrates as rough as possible. We have shown that, for a 
given solid fraction and above a roughness threshold the value of the energy barrier stabilizing 
the fakir states remains constant. Thus, putting efforts to design ultra rough surface seems 
ineffective. Conversely, reducing the size of the elementary pattern of the surface would arbi- 
trarily increases the resistance against impalement [21]. We precise that this conclusion relies 
on the quantitative agreement between our impact and quasistatic experiments. However, for 
impacts on substrates patterned at ultra-small scales, the confinement of the fluid flows would 
enhance the magnitude of the hydrodynamic forces hindering impalement. The identification 
of this characteristic scale at which they would overcome capillary effects remains an open 
question. 
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